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Abstr act. Relativistic thermodynamics is constructed from the point of view 
of special relativistic hydrodynamics. A relativistic four-current for heat and a 
general treatment of thermal equilibrium between moving bodies is presented. 
The different temperature transformation formulas of Planck and Einstein, 
Ott, Landsberg and Doppler appear upon particular assumptions about inter- 
nal heat current. 



1. Introduction 

Considering the temperature of moving bodies, the easier question is to answer, 
what is the apparent spectral temperature. In this case a spectral parameter is 
transformed if the thcrmalized source is moving with respect to the observer (de- 
tector system), and the transformation rule can be derived from that of the energy 
and momentum in the co-moving system. 

This has been known from the beginnings of the theory of special relativity and 
never has been seriously challenged. An essentially tougher problem is to under- 
stand the relativistic thermalization: what is the intensive parameter governing 
the state with energy exchange equilibrium between two, relatively moving bodies 
in the framework of special relativity. In particular how this general temperature 
should transform and how does it depend on the speed of the motion. Here several 
answers has been historically offered, practically including all possibilities. 

Planck and Einstein concluded that moving bodies are cooler by a Lorentz factor 
[U [21 [3] , first Blanusa then Ott has challenged this opinion [H [5] by stating that 
on the contrary, such bodies are hotter by a Lorentz factor. During later disputes 
several authors supported one or the other view (see e.g. [6j [7j QUI HJJ H21 
ITS"] and the references therein) and also some new opinions emerged. Landsberg 
argued for unchanged values of the temperature Q31 [15] . Other authors observed 
that for a thermometer in equilibrium with black body radiation the temperature 
transformation is related to the Doppler formula [16] [TT1 [TBI El EO] > therefore the 
measured temperature seems to depend on the physical state of the thermometer. 
This problem is circumvented by the suggestion that thermal equilibrium would 
have a meaning only in case of equal velocities [21] [22] [23] . 

Behind these different conclusions there are, in our opinion, different views about 
the energy transfer and mechanical work, and the identification of the heat [12] . In 
a simplifying manner the assumptions and views about the Lorentz transformation 
properties of internal energy, work, heat, and entropy influence such properties 
and the very definition of the absolute temperature. Coming to the era of fast 
computers, a renewed interest emerged in such questions by modelling stochastic 
phenomena at relativistic energy exchanges and relative speeds [24J [25] [26]. In 
particular, dissipative hydrodynamics applied to high energy heavy ion collisions 
requires the proper identification of temperature and entropy (27] [28j [29l [30l EJJ 
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132] 133] 134] [35] . In this letter we show that our approach to replace the Israel- 
Stewart theory of dissipative hydrodynamics, proposed earlier [36] [32 [38l EH], is 
related to the problem of thermalization of relatively moving bodies with relativistic 
velocities and our suggestion is compatible with the foundations of thermodynamics 
and guarantees causal heat propagation. 

00 

By doing so we encounter the following questions in our analysis: 

(1) What moves (or flows)? Total energy and momentum do flow correlated, 
but further conserved charges (baryon number, electric charge, etc.) may 
flow differently. In relativistic systems one has to deal with the possibility 
that the velocity field is not fixed to either current, not being restricted to 
the Landau-Lifshitz [40] or Eckart [41] frames. 

(2) What is a body? We exploit, how do integrals over extended volumes relate 
to the local theory of hydrodynamics, and what is a good local definition 
for volume change in relativistic fluids. In close relation to this, we suggest 
a four- vector generalization to the concept of heat. 

(3) What is a proper equation of state? Here the functional dependency be- 
tween entropy and the relativistic internal energy is fixed to a particular 
form. 

(4) What is the proper transformation of the temperature? As we have men- 
tioned above prominent physicists expressed divergent opinions on this in 
the past. This problem is intimately related to that of thermal equilibrium 
and to the proper description of internal energy. 



2. Hydro- and thermodynamics 

In this letter we concentrate on the energy-momentum density of a one-component 
fluid, but the results can be generalized considering conserved currents in multi- 
component systems easily. The energy-momentum tensor can be split into compo- 
nents aligned to the fiducial four- velocity field, u a (x) , and orthogonal ones: 

(1) T ab = eu a u b + u a q b + q a u b + P ab 

with u a q a = and u a P ab = P ab u b = 0. When considering complex systems, like a 
quark-gluon plasma, the velocity field can be aligned only with one of the conserved 
currents, unless several currents are parallel (i.e. different conserved charges are 
fixed to the same carriers) . In our present treatment the velocity field is general. 

Relativistic thermodynamics is obtained by integrating the local energy-momentum 
conservation on a suitably defined extended and homogeneous thermodynamic 
body. Therefore in the balance of energy-momentum we separate the terms per- 
pendicular and parallel to the velocity field as 

d b T ab = ^-(eu a +q a ) + (eu a + q a )d b u b 
or 

+ p(±u a + u a d b u b )-^(u a q b + Il ab ) 

(XT CLT 

(2) - V> + V b (u a q b + U ab ) = a . 

From now on the proper time derivative is denoted by a dot / = df /dr = u a d a f 
for an arbitrary function f(x). V a = d a ~u a u c d c denotes a derivative perpendicular 
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to the velocity field and we also split the pressure tensor into a hydrostatic part 
and a rest: P ab = p(u a u b - g ab ) + Il ab . 

Let us now assume, that u a is smooth and we may give a connected smooth 
surface H that is initially perpendicular to the velocity field and has a smooth (two 
-dimensional) boundary. As a further simplification we will assume that the velocity 
field is not accelerating u a = a , therefore d a u a — V Q u a and the hypersurface 
remains perpendicular to the four-velocity field. Hence the propagation of the 
surface can be characterized by the proper time r of any of its wordlines. We 
refer to this hypersurface - a three dimensional spacelike set related to our fluid - 
as a thermodynamic body. Considering homogeneous bodies we set V a e = and 
V a p = 0. It is important that the velocity field itself is not homogeneous, V a u b ^ 0. 
Now integration of ([2]) on H(r) results in 

{eu a + q a + (eu a + q a )d b u b + pu a d b u b )dV = 

H{r) 

(3) / V b (u a q b +n ab )dV. 

Jh(t) 

With the above conditions wa apply the transport theorem of Reynolds to the l.h.s. 
of eq.© and the Gauss-Ostrogradsky theorem to the r.h.s. of eq.Q and obtain 

(4) Eu a + G a + P u a V = I (u a q b +n ab ) dA b = 5Q a . 

dH{r) 

Here u a — f„ u a dV/V is the average velocity field inside H , E — eV is the total 
energy, G a — J H q a dV, and dA b is the two-form surface measure circumventing the 
homogeneous body in the region H(t). The two-dimensional surface integral term 
is the physical energy and momentum leak (dissipation rate) from the body under 
study, we denote it by SQ a . This is a four- vector generalization of the concept of 
heat. It describes both energy and momentum transfers to or from the homogeneous 
body. 

The derivation of the temperature in thermodynamics is related to the maximum 
of the total entropy of a system (under various constraints). This way its reciprocal, 
1 /T is an integrating factor to the heat in order to obtain a total differential of the 
entropy [42 , 43. . Here we follow the same strategy considering a vectorial integrating 
factor A a : 

(5) 5Q a = E a + pu a V = A a S + S a 

with E a — Eu a + G a the energy-momentum vector of the body, u a = and 
S a orthogonal to A a . The decisive point is, that - according to the above - the 
entropy of the homogeneous body is a function of the energy-momentum vector and 
the volume: S = S(E a 7 V). Multiplying eq.© by A a dT/(A b A b ) and utilizing that 
du a — we obtain 

The connection to classical thermodynamics is best established by the temperature 
definition 

1 A a u a 



(7) 



T ' A b A b 
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The intensive parameter associated to the change of the four-vector G a is denoted 

by 

g a A a 

With these notations we arrive at the following form of the Gibbs relation: 

(9) TdS = dE + g a dG a + pdV. 

Due to the definitions eq. (|7I8[> g a u a = 1. Hence the Gibbs relation can be written 
in the alternative form 

(10) TdS — g a dE a + pdV, 

suggesting that the traditional change of the energy, dE, has to be generalized to 
the change of the total energy-momentum four- vector, dE a = d(Eu a + G a ). 

For the well-known Jiittner distribution [33] g a = u a . This equality has been 
postulated among others in the classical theory of Israel and Stewart [35]. Then 
(fTU]) reduces to 

(11) TdS = dE + pdV, 

where E = u a E a . In this case the internal energy can be interpreted as E, but 
its total differential contributes to the Gibbs relation. In the general case g a ^ u a 
- considered below - there remains a term, related to momentum transfer. It is 
reasonable to assume that the new intensive variable is timelike: g a g a > 0. Then 
we introduce 

(12) w a :=g a -u a . 

Now ||u> a || 2 = —w a w a = — g a g a + l<l follows. The spacelike four-vector w a has the 
physical dimension of velocity. Due to — 1 < w a w a < and u a w a = its general 
form is given by w a = (7u|w|,7w). In this case |w| 2 < 1. We interpret w as the 
velocity of the internal energy current. 

Here some important physical questions arise: is it only a single or several differ- 
ential terms describing the change of energy and momentum? When two, relatively 
moving bodies come into thermal contact what can be exchanged among them in 
the evolution towards the equilibrium? 



3. TWO BODIES IN EQUILIBRIUM 

Let us now consider two different bodies with different average velocities and 
energy currents. When all components of E a and the total volume are kept con- 
stant independently, i.e. dE? + dE% = and dV x + dV 2 = while dS(Ef,Vi) + 
dS(E%, V2) = in the entropy maximum, then from (|10p we obtain the conditions 

(13) gl = gi pi _ n 

1 ' Ti T 2 ' li T 2 ' 

This, in general, does not mean the equality of temperatures. 

In order to simplify the discussion we restrict ourselves to one-dimensional mo- 
tions and consider u a = (7,71") with 7 = 1/yl — v 2 Lorentz-factors. The energy 
current velocity is given by w a = (jvw^w) and q a = (7(1 + vw),^(v + w)). Here 
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w describes the speed of internal energy current. The thermal equilibrium condition 
(f]~3|) hence requires 



7i(l + wiWi) _ 72(1 + ^2^2) 

(14) 



Tt T 2 
7i(«i+wi) 72(^2+^2) 



Ti T 2 

The ratio of these two equations reveals that in equilibrium the composite relativis- 
tic velocities are equal, 

, 15 v vi + Wi _ V2 + W2 

1 + V\W\ 1 + V 2 W 2 ' 

and the difference of their squares leads to 



(16) • 1 - • 



Ti T 2 

The equality of some other velocities were investigated by several authors [211 E21 

cum]. 

One realizes that in the thermal equilibrium condition four velocities are involved 
for a general observer: v\, v 2 , wi and w 2 . By a Lorentz transformation only one of 
them can be eliminated. The remaining three (relative) velocities reflect physical 
conditions in the system. According to eq. (|15|) 



v + w 2 

17 wi = — 

1 + vw 2 



with v = (v 2 — v%) / (1 — V1V2) relative velocity. The associated factor, yT — w± can 
be expressed and the temperatures satisfy 



(18) Ti = T 2 



1 + vw 2 

This includes the general Doppler formula [TBI IT71 U51 \W[ \M ITT] . 

It is enlightening to investigate this formula with different assumptions about the 
energy current speed in the observed body, w 2 . The induced energy current speed 
in an ideal thermometer, wi and the temperature it shows, Ti, are now determined 
by eas. (fl~7)) and ([PS)) . Figure [JJ plots temperature ratios T\jT 2 for a body closing 
with v = —0.6 as a function of the energy current speed, w 2 - 

(1) w 2 = 0: the current stands in the observed body. In this case w\ = 
v, the measured energy current speed is that of the moving body, and 
Ti = T 2 ^fY^ 1 < T 2 , the moving body appears cooler by a Lorentz factor 
CD El 03 (see Fig © . 

(2) wi — 0: the current stands in the thermometer. In this case we must 
have u>2 = —v and T\ = T 2 j\J\ — v 2 > T 2 , the moving body appears hotter 
CHE1CDCL2] (see Fig©. 

(3) wi + w 2 = 0: the current is standing in the total system of moving body 
and thermometer, the individual contributions exactly compensate each 
other. This is achieved by a special value of the energy current velocities, 
w 2 = —w, w\ = w with w — (1 — vl — v 2 ) /v. In this case even the apparent 
temperatures are equal, T\ = T 2 CHI CCS] (see Fig [4]). 
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(4) W2 = 1: a radiating body (e.g. a photon gas) is moving. In this case 

wi = 1, and one obtains Tj = 22^/ ■jfj^j. It means that T\ < T2 for v > 0, 

a Doppler red shifted temperature is measured for an aparting body (see 
Fig EJ - quite common for astronomical objects - and T\ > T2 for v < 0, a 
Doppler blue shifted temperature appears for closing bodies - more common 
in high energy accelerator experiments. 

On figures ©-© we fix the reference frame to the thermometer, therefore 
uj = (1,0) (the vertical axis is time). The energy current velocity four- vectors 
are perpendicular to the corresponding four-velocities, therefore they are on lines 
symmetrical to the light cones. The four-velocity vectors end on the timelike hy- 
perbolas and the spacelike energy current velocities end inside the spacelike hyper- 
bolas. The temperature ratios are determined by the magnitudes of the g a -s as 
Ti/T 2 = \\gn/m (see 08]. 

4. LORENTZ SCALAR TEMPERATURE 

According to the classical ansatz g a = u a , the total entropy has to depend on 
the total energy E = u a E a . Then the equilibrium conditions (TTTl) and (fl8|) result 
in zero relative velocity v = and the temperatures are equal T± — T 2 . 

However, thermodynamic and generic stability considerations are favoring an 
other Lorentz-scalar combination \\E\\ = y/E a E a [36] |3Zj [38l [39]. Denoting the 
partial derivative of entropy 5(||.E||, V) with respect to its first argument by 1/0 = 
, one re- writes the total differential, 

(19) dS = ±d\\E\\ + V -dV = ^dE a + ^dV, 

and comparing to the general Gibbs relation (jlOl) one obtains the correspondence 

g a 1 E a P P 



(20) 



T e\\E\v T 



It follows that the length of the intensive four- vector, g a , is the ratio of the tradi- 
tional (energy associated) and scalar (energy-momentum four- vector length associ- 
ated) temperatures: ^Jg a g a = T/6. On the other hand its projection to the average 
velocity reveals the value in the comoving system: 

This equation relates the energy-momentum-associated scalar temperature, 9 to 
the energy-associated one, T. As a consequence we obtain 

(22) g° = 41, = Ea -^ Eb ^. 

The later formula clearly interpret w a as the quotient of the comoving, average 
velocity related energy current (momentum) and energy of the thermodynamic 
body, that is the energy current velocity. 

Finally we remark, that in the simple two dimensional particular case we obtain 
that 8 = T I \J\ — w 2 . Therefore the equilibrium condition (fl6]) gives equal scalar 
temperatures: Q\ =62- This is a stronger reflection of Landsberg's view, and his 
physical arguments in [141 115] than the assumption of zero total energy current 
velocity. 
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Figure 1 . Ratio of the temperatures of the observed body in its 
rest frame, T2 to that shown by an ideal thermometer, T\ as a 
function of the the speed of the heat current in the body, W2 while 
approaching with the relative velocity v = —0.6. 




Figure 2. The space-time figure for the Planck-Einstein rule of 
two thermodynamic bodies in equilibrium. There is no energy 
current in the observed body (w^ = 0), therefore the four- 
velocity (solid arrow) is parallel to the vectors (g^g^)- The ratio 
of the temperatures is Ti/T 2 < 1. 

5. Summary 

We investigated the possible derivation of basic thermodynamical laws for ho- 
mogeneous bodies from relativistic hydrodynamics. The dependence of entropy 
on internal energy is replaced by a dependence on the energy-momentum four- 
vector, E a . As a novelty a relativistic heat four- vector has been formulated. For 
the traditional, energy exchange related temperature, T, a universal transforma- 
tion formula is obtained. For a general observer four velocities are involved in the 
equilibrium condition of two thermodynamic bodies in equilibrium. One of them 
can be eliminated by choosing the observing frame, the physical relation depends 
only on the relative velocity. Another condition connects the internal heat currents 
in the bodies in thermal contact. So there remains two velocity like parameters 
to describe thermal equilibrium: the energy current speed (the velocity related to 
the integrated internal heat current density) in one of the bodies and their relative 
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Figure 3. The space-time figure for the Blanusa-Ott rule. The 
energy current stands in the thermometer (wf = 0), therefore the 
Ui four- velocity (solid arrow) is parallel to the vectors (gf, g$). The 
ratio of the temperatures is T\jTi > 1. 




Figure 4. The space-time figure for the Landsberg rule. There 
is no energy current in the composed system, therefore the four- 
vectors (<7°, g% , dotted arrows) are equal. The ratio of tempera- 
tures is Ti/T 2 = 1. Here w x = -0.33, w 2 = 0.33. 



velocity. The traditional temperature transformation formulas belong to corre- 
sponding particular choices on the energy current speeds. This can be the reason 
that no agreement could be achieved historically. For most common cases there is 
no heat current in the observed body but it flows in the thermometer. This leads 
to the Planck-Einstein transformation formula. 

The closer relation to dissipative hydrodynamics favors a particular dependence 
of entropy on energy-momentum and leads to a Lorentz scalar temperature. 

Our approach is covariant, and the compatibility to hydrodynamics clarifies that 
the Planck-Ott imbroglio is not a problem of synchronization as it was supposed in 
[121I2S]. It makes possible to interpret the classical paradoxical results of Planck and 
Einstein, Ott, Landsberg and Doppler in a unified treatment. Our investigations 
reveal that despite of the apparent paradoxes related to Lorentz transformations, 
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Figure 5. The space-time figure for the Doppler red shift rule 
of two thermodynamic bodies in equilibrium. The energy cur- 
rent speed (dashed arrows) in the observed body is that of the 
light W2 = 1, therefore wi = 1, and the four- vectors of energy- 
momentum intensives g%, dotted arrows) are light-like. 

there is a covariant relativistic thermodynamics with proper absolute temperature 
in full agreement with relativistic hydrodynamics. 
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